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and Redshift Data. II. Results from an Expanded Sample 
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ABSTRACT 



This is the second in a series of papers in which we compare Tully-Fisher (TF) data 
from the Mark III Catalog with predicted peculiar velocities based on the IRAS galaxy 
redshift survey and gravitational instability theory, using a rigorous maximum likelihood 
method called VELMOD. In Paper I (Willick et al. 1997b), we applied the method to a 
cz LG < 3000 kms" 1 , 838-galaxy TF sample and found ft = 0.49 ± 0.07, where ft = fi°- 6 /&i 
and bi is the linear biasing parameter for IRAS galaxies. In this paper we increase the 
redshift limit to czlg — 7500 kms -1 , thereby enlarging the sample to 1876 galaxies. The 
expanded sample now includes the W91PP and CF subsamples of the Mark III catalog, in 
addition to the A82 and MAT subsamples already considered in Paper I. 

We implement VELMOD using both the forward and inverse forms of the TF relation, 
and allow for a more general form of the quadrupole velocity residual detected in Paper I. 
We find ft = 0.50 ±0.04 (la error) at 300 kms -1 smoothing of the IRAS-predicted velocity 
field. The fit residuals are spatially incoherent for ft = 0.5, indicating that the IRAS plus 
quadrupole velocity field is a good fit to the TF data. If we eliminate the quadrupole we 
obtain a worse fit, but a similar value for ft of 0.54 ± 0.04. Changing the IRAS smoothing 
scale to 500 kms -1 has almost no effect on the best ft. We find evidence for a density- 
dependence of the small-scale velocity dispersion, a v (5 g ) ~ (100 + 35 5 g ) kms -1 . 

We confirm our Paper I result that the TF relations for the A82 and MAT samples 
found by VELMOD are consistent with those that went into the published Mark III catalog. 
However, the VELMOD TF calibrations for the W91PP and CF samples place objects ~ 8% 
closer than their Mark III catalog distances, which has an important effect on the inferred 
large-scale flow field at 4000-6000 kms -1 . With this recalibration, the IRAS and Mark III 
velocity fields are consistent with one another at all radii. 



1. Introduction 

In recent years, a number of groups have compared the peculiar velocity and/or density fields derived 
from distance indicator data with the corresponding fields obtained from redshift survey data (Kaiser et 
al. 1991; Dekel et al. 1993; Hudson 1994; Roth 1994; Hudson et al. 1995; Schlegel 1995; Davis, Nusser, 
& Willick 1996; Willick et al. 1997b, hereafter Paper I; da Costa et al. 1997; Riess et al. 1997; Sigad 
et al. 1998). The principal goals of these comparisons are to test the gravitational instability (GI) 



picture for the growth of large-scale structure and to measure the parameter ft = Q os /b, where Q is 
the present value of the cosmological density parameter, and b is the "biasing parameter" (see below). 
A longer-range goal is to measure Q itself, by combining the ftmeasurement with other measurements 
that constrain a combination of Q and b. 

Measurement of ft is based on the relationship between the peculiar velocity and density fields pre- 
dicted by GI for the linear regime (Peebles 1980): 



In equation (JI|), the galaxy number density fluctuation field 5 g is assumed to be related to the underlying 
mass density fluctuation field S by the simple linear biasing model S g = bS. Taking the divergence of 
both sides of equation ([!]) yields: 

V-v p = -ft6 g . (2) 

In both equations, distances are assumed to be measured in units of the Hubble velocity (i.e. Hq = 1). 

To estimate the quantity ft via equation ([]]), one measures S g from redshift survey data, and then 
predicts v p (r) for a sample of galaxies with redshift-independent distances^ and thus estimated peculiar 
velocities. One then asks, for what value of ft does the velocity field prediction best fit the TF data? 
This approach is known as the "velocity-velocity" (v-v) comparison. 

Alternatively, one can do a "density-density" (d-d) comparison, using equation (0). In this case the 
crucial input from the redshift survey is not the predicted velocity field v p (r), but instead the directly 
observed density field 5 g . However, the TF data must now be converted into a three-dimensional velocity 
field, whose divergence is then taken to yield an effective mass density field —V • v p . Comparison of 5 g 
and — V ■ v p , via equation (Q), then yields ft. 

In the v-v comparison the redshift survey data is manipulated to yield predicted peculiar velocities 
(see, e.g., Yahil et al. 1991). The way this predicted velocity field changes with ft is what provides the 
v-v comparison with its discriminatory power. In the d-d comparison it is the numerical processing of 
the TF data which is more important for /^-determination. This is done using the POTENT method 
and its variants (Bertschinger & Dekel 1989; Dekel, Bertschinger, & Faber 1990; Dekel 1994, 1997; da 
Costa et al. 1996) which invoke the assumption of potential flow in order to convert the radial TF data 
into a 3-dimensional velocity field, and thus into an effective mass density field. 

The redshift survey most often used in recent v-v and d-d comparisons, and the one we use in this 
paper, is the 1.2 Jy IRAS redshift survey (Fisher et al. 1995), which covers nearly the full sky and is 
only weakly affected by dust extinction and related effects at low Galactic latitude. Hereafter, we write 
bj to denote the IRAS biasing parameter, and ft = Q os /bi. 

The published results for ftj appear to bifurcate according to whether the d-d or the v-v comparison 
is used. The former has been implemented using the POTENT method by Dekel et al. (1993), Hudson 
et al. (1995), and Sigad et al. (1998; hereafter POTIRAS) to obtain ft = 1.29 ± 0.30, ft ~ 1.0 ± 0.13, 
and ft = 0.89 ± 0.12, respectively^ (the error bars are la). In the first of these studies, POTENT was 



1 Hereafter we will assume for definitcncss that the redshift-independent distances have been derived from the Tully- 
Fisher (1977; TF) relation. However, the method applies to any distance indicator relation. 

2 The Hudson et al. result has been converted from the measured (3 op t assuming that b op t/bi = 1.3 (Strauss et al. 1992). 
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applied to the redshift-independent distances in the Mark II Catalog (Burstein 1989), while in the latter 
two it was applied to those in the Mark III Catalog (Willick et al. 1997a). These relatively high values 
of j3[ have often been cited (assuming the that bj is not much different from unity) as evidence for an 
Q = 1 universe. In contrast, the v-v approach has typically produced lower values of Pi, which (again 
assuming that bj « 1) point to a low-density (Q ~ 0.2-0.5) universe. Davis, Nusser, & Willick (1996) 
and da Costa et al. (1997) each found Pi = 0.6±0.15 by applying the inverse Tully-Fisher (ITF) method 
of Nusser & Davis (1995), in the former case to the Mark III catalog and in the latter case to the SFI 
sample of Giovanelli et al. (1997). Riess et al. (1997) also used the ITF method for distances obtained 
from Type la supernovae, finding f3j = 0.40 ± 0.15. Roth (1994) and Schlegel (1995) used v-v analyses 
of smaller TF samples to obtain (3j = 0.6 and Pi = 0.39 respectively. Shaya, Peebles, & Tully (1995) 
find (3j = 0.45 ± 0.15 from their v-v analysis of nearby TF dataQ. Finally, we found (3j = 0.49 ± 0.07 
in Paper I by applying the VELMOD method (§ 2) to a subset of the Mark III Catalog restricted to 
cz LG < 3000 km s _1 . See Strauss & Willick (1995, hereafter SW) for a review of these and other methods 
for measuring /3. 

In this paper we will again apply VELMOD, now to an expanded sample that includes all Mark III 
Catalog field spirals out to czlg — 7500 kms -1 . This larger sample will lead to tighter constraints on 
Pi than obtained in Paper I, although our results will be fully consistent. The outline of this paper is 
as follows. In § ^| we review the VELMOD method. In § |3| we describe the selection of our expanded 
sample. In § f|, we discuss the motivation behind and implementation of a more general form of the 
quadrupole velocity residual introduced in Paper I. In § [|, we present the main results of the maximum 
likelihood analysis. In § [], we compare the VELMOD TF calibrations to those in Mark III. In § 0, we 
quantify the goodness of fit of our model to the data. Finally, in § |] we summarize our main conclusions. 
In the Appendix, we describe an analytic approximation to computing the VELMOD likelihoods. The 
formulae presented there are not limited to VELMOD and are generally useful in velocity field analyses. 



VELMOD is a maximum likelihood method for comparing TF data to predicted peculiar velocity 
fields. The method was described in some detail in Paper I, §2, and we give only a brief overview 
here. The TF data for each galaxy consist of its direction (/,&), its redshift cz measured in the Local 
Group (LG) frame, its apparent magnitude m, and its velocity width parameter i] = \og(Av) — 2.5. 
The velocity field model gives the relationship between redshift and distance (r) along any given line 
of sight, albeit with some finite scatter, <j v , due to inaccuracies of the model and small-scale velocity 
"noise." We assume that there exists a forward [M(??)] and an inverse [?7 (M)] TF relation for each 
sample, such that m, 77, and r are related as follows: 



2. Method of Analysis 
2.1. The VELMOD Approach 



m = M(r]) + 5 logr 



A — b i] + 5 log r , with rms dispersion <7tf 



(3) 



(forward relation), or 



i] = rf{m — 5 logr) 



e (m — 5 logr — D) , with rms dispersion a v 



(4) 



3 We have again converted their f3 opt into an equivalent /3/ assuming b opt /bi = 1.3. 
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(inverse relation). We refer to A, b, and o"tf (D, e, and a v ) as the zero point, slope, and scatter of the 
forward (inverse) TF relation, or simply as the TF parameters. 

For each object in the TF sample, P(m\r], cz) — the probability that a galaxy of redshift cz and velocity 
width parameter rj will have apparent magnitude m — is evaluated when the forward TF relation is used; 
see equation (|A1|). For the inverse TF relation, it is P(r]\m, cz) that is evaluated, equation (|A^). These 
single-object probabilities depend on a number of parameters: 

1. The three TF parameters for each distinct subsample (i.e., A82, MAT, W91PP, and CF). In 
§ |5.1] , we explore the addition of a fourth TF parameter describing the change in the scatter with 
luminosity. 

2. j3i, which determines the IRAS-predicted peculiar velocity. 



3. The small-scale velocity dispersion o~ v . In § |5.4j , we include an additional parameter f$ describing 
the density dependence of a v . 

4. A cutoff scale, Rq, for the velocity quadrupole (§ §). 

5. A LG velocity vector wlg, required because small errors in the prediction of the LG velocity 
propagate to all other peculiar velocity predictions (cf. Paper I, § 2.2.3). As wlg is primarily 
determined by nearby galaxies, in this paper we simply fix it to its Paper I value. 

The single-object probabilities are multiplied together, yielding an overall probability P for the entire 
TF sample. The value of (3i for which P is maximized is the maximum likelihood value of (3j. In practice, 
rather than maximizing P we minimize C = —2 In P. (In § 5, we will write £f or w and C llw to distinguish 
forward from inverse likelihoods.) A single VELMOD run consists of minimizing C, at each of 10 values 
of /3j, 0.1, 0.2, . . . , 1.0, by continuously varying the TF parameters of each sample?]. A cubic fit to the 
£(/?/) points then yields the maximum-likelihood value of /?/. Tests with mock catalogs, discussed in 
Paper I, demonstrated that this maximum likelihood value of (3i is an unbiased estimator of the true 
value when the IRAS peculiar velocities are predicted using a 300 kms -1 Gaussian smoothing scale and 
a Wiener filter. The tests also showed that rigorous 1 a errors in f3i are given by noting the values at 
which L differs by one unit from its minimum value, as obtained from the cubic fit. 

Because the TF parameters for each sample are determined via maximum likelihood, a priori TF 
calibrations are not required for VELMOD. Indeed, each value of /?/ is given the fairest possible chance 
to fit the data by finding the TF parameters most in accord with the velocity field it produces. These 
TF parameters are not constrained to be similar to those used to produce the Mark III catalog distances 
(we discuss this issue further in § ^j). Furthermore, while the TF scatter is treated as a free parameter, 
we emphasize that maximizing likelihood is not equivalent to minimizing scatter (cf. Paper I, § 3.4). In 
general, the minima of £ and of <7tf (or a v ) for a given subsample do not precisely coincide. 



4 We hold the velocity parameters a v and Rq fixed in any given VELMOD run, but carry out a series of runs in which 
they take on a range of discrete values, and in this way determine their maximum likelihood values; cf. § 5.3 and 5.4. The 
only velocity parameter treated as continuously variable is fg. 
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2.2. Implementation of Inverse VELMOD 

Because selection effects on the forward TF relation are strong (Willick 1994), the sample selection 
function must be properly modeled in forward VELMOD in order to obtain unbiased results. However, 
as selection depends weakly on velocity width, errors in modeling the selection function will have little 
effect on inverse VELMOD or comparable analyses. For this reason, inverse TF methods have been 
favored by many workers (e.g., Schechter 1980; Aaronson et al. 1982b; Tully 1988; Nusser & Davis 1995; 
Shaya et al. 1995; da Costa et al. 1997; cf. SW for a discussion). On the other hand, inverse, but not 
forward, VELMOD depends on the galaxy luminosity function $(M) (cf. Paper I, § 2), which is not 
easy to quantify, given the fact that each sample uses its own photometric system. 

Because $(M) appears in the integrals in both the numerator and denominator of the expression 
for P(j)\rn, cz) (equation ||A2|| ), it is not crucial to model it perfectly. We determine $(M) for each 
sample as follows. As 77 is defined in essentially the same way for each sample, we assume that there 
is a universal ^-distribution function, 4>(t]), which we take to be a Gaussian of mean 770 = —0.05 and 
dispersion = 0.15. This distribution function matches well what is seen in the Mark III TF samples 
above the cutoff r/ m i n ~ —0.4 imposed by magnitude and diameter limit effects. We then calculate $(M) 
using the relationship between <f)(r]) and $(M) given by the TF relation itself: 



$(M) 



77°(M) 



dM 



e<p\r] {M)] , (5) 



where rj (M) is the inverse TF relation and e is its slope (cf. equation |4[]). 

The luminosity function obtained from equation @ is, as required, different for each sample, because 
each sample has its own TF parameters. The differences reflect bandpass effects and differing approaches 
to extinction/inclination corrections for each of the individual Mark III TF samples (Willick et al. 1997a). 
Ultimately, we will test the suitability of this approximation by comparing the results of the forward 
and inverse VELMOD calculations. To the extent they agree, we can be confident that our imperfect 
modeling of the selection and luminosity functions do not bias the results. 



2.3. An Analytic Approximation to the VELMOD Likelihoods 

A drawback of the original VELMOD algorithm was its repeated evaluation of the numerical integrals 
in terms of which the single-object likelihoods are defined. These integrations are crucial in triple- 
valued or flat zones in the redshift-distance relation (cf. Paper I, § 2.2.2). However, away from such 
regions, and at distances much larger than a v , maximizing the VELMOD likelihood is very similar 
to minimizing differences between TF distances and those inferred from the velocity field model (the 
"Method II" approach to velocity analysis; SW, § 6.4.1). This suggests that we can find an accurate 
analytic approximation to the exact VELMOD likelihoods for many galaxies. Equation (15) of Paper 
I is an approximation for the forward likelihood in the simple case when selection effects are neglected 
and a constant density field is assumed. We have since generalized this result to all relevant cases and 
applied it in our calculations, thereby reducing the run time of the code by a factor of ~ 4. The details 
are complex and are given in the Appendix; we discuss the salient features here. 

For each TF galaxy, the velocity field model yields a "crossing-point" distance w, defined implicitly 
by w + u(w) = cz, where u(r) is the radial component of the predicted peculiar velocity along the line 
of sight. Similarly, the TF relation defines a distance d implicitly by 51og<i = m — M(t]) (forward) or 
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r] (m — 51og<i) = 7] (inverse). Our main result is that the forward and inverse single-object likelihoods 
are well approximated as normal distributions in \n(d/w), but with the mean value of ln<i offset from 
\nw by an amount proportional to A^, where A„ = a v /[w(l + u')}, and u' = [du/dr] r=w . At large 
distances, w ^> a v , and for u' ^,0, A^ is very small and our approximation becomes more accurate. In 
practice we have found that the velocity field is cold (a v ~ 130 kms -1 ; § p.4[ ), and A^ is usually small 
even at distances as low as 1000 kms -1 . 

We have checked the analytic approximation against the full numerical integration in a number of 
regimes. It fails in regions where the velocity field changes sufficiently rapidly in the vicinity of the 
crossing point. In practice, this invalidated the approximation in a 30° cone around the Virgo cluster, 
for w < 2700 kms -1 . More generally, we found that it is inaccurate for A„ > 0.2 regardless of the nature 
of the velocity field. The likelihoods for all such objects were always computed using the full numerical 
integration. However, for the remaining ~ 75-80% of sample galaxies, we found that the approximation is 
remarkably accurate. The rms difference between the exact and approximate likelihoods for such objects 
is 0.015 in C. This accuracy is sufficient to minimize L at a given f3i by varying the TF parameters and 
relevant velocity parameters other than /3j. Once this minimum is found, we re-evaluate C using the 
exact numerical probabilities for all objects. The final maximum-likelihood value of (3i is derived from 
these exact values of C. However, this maximum-likelihood value always differed by < 0.01 from one 
obtained from the approximations. Thus, we are confident that our use of the approximate likelihoods 
in the parameter variation procedure has not affected our maximum likelihood results for (3i. 

3. Selection of the Expanded Sample 

In Paper I we limited our analysis to the local (czlg < 3000 kms -1 ) volume. This constrained us to 
use only two TF subsamples of the Mark III Catalog: the Aaronson et al. (1982a; A82) and Mathewson 
et al. (1992; MAT) data sets. The former is a 1.6/xm (H band) photometry, 21 cm velocity width data 
set; the latter consists of / band CCD magnitudes and a mixture of 21 cm and optical velocity widths 
(cf. Willick et al. 1997a for further details). The remaining Mark III TF samples contain too few galaxies 
within 3000 kms -1 to have made their inclusion worthwhile. 

Here we increase our redshift limit to c^lg — 7500 kms -1 . However, because the A82 sample itself 
is badly incomplete beyond 3000 kms -1 (Willick et al. 1996), we continue to use the same 300 galaxy, 
c^lg < 3000 kms -1 A82 subsample used in Paper I. Our MAT TF sample, in contrast, has grown from 
538 galaxies in Paper I to 1159 galaxies for the present analysis. No changes were made in the way we 
select the MAT galaxies. Specifically, we continue to apply a (photographic) diameter limit of 1.6 arcmin 
and to require that log(a/6) > 0.1, where a/b is the major to minor axis ratio. The last requirement 
excludes objects that are too face-on and which thus have large velocity width uncertainties. 

With the higher redshift limit we now include the two other TF field samples in the Mark III catalog, 
the Willick (1991; W91PP) Perseus-Pisces sample and the Courteau-Faber (Courteau 1992, 1996; CF) 
Northern sky sample. Both of these data sets consist of R band CCD magnitudes. W91PP uses the 
21 cm velocity widths of Giovanelli et al. (1985, 1986) and Giovanelli & Haynes (1989), while CF 
uses optical velocity widths (Courteau 1997). W91PP and CF were originally designed to have uniform 
photometric and velocity width properties. The mutual consistency of the photometry for these samples 
has indeed been verified (Willick 1991; Courteau 1992, 1996). However, their velocity widths have not 
been shown to be consistent, and Willick et al. (1996, 1997a) found different TF calibrations for the two 
samples, as we will here. 
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TABLE 1 

Mark III Subsamples Used in the VELMOD Analysis 



Sample 


Redshift limit 


Mag/Diam Limit 


N 


notes 


A82 


3000 kms" 1 


14.0 mag 


300 


a,b 


MAT 


7500 kms" 1 


1.6 arcmin 


1159 


c 


W91PP 


7500 kms" 1 


1.0 arcmin 


247 


d 


CF 


7500 kms" 1 


1.5 arcmin 


170 


e 


Total 


1876 



Table 1: Notes: (a) The A82 sample used in this paper is identical to that used in Paper I; it is very 
incomplete beyond 3000 kms -1 and hence was not extended to a higher redshift limit, (b) The A82 
magnitude limit applies to the RC3 catalog blue magnitudes; see Willick et al. (1996) for further details, 
(c) The MAT diameter limit applies to the ESO catalog blue photographic diameters; see Willick et al. 
for further details, (d) The formal diameter limit used in calculating the W91PP selection function was 
1.15 arcmin; the diameters in question are UGC blue photographic diameters, (e) The formal diameter 
limit used in calculating the CF selection function was 1.6 arcmin; UGC blue photographic diameters 
are again used. 



The selection criteria for W91PP and CF are not known as rigorously as might be hoped. Both 
samples are selected to the limit of the UGC catalog — nominally, therefore, to a photographic diameter 
limit of 1.0 arcmin. However, the UGC catalog is known to become increasingly incomplete below about 
1.5 arcmin (Hudson & Lynden-Bell 1991). This problem was studied by Willick et al. (1996), who found 
that consistency of W91PP group distance moduli as measured by the forward and (essentially selection- 
bias free) inverse forms of the TF relation was achieved with a diameter limit of 1.15 arcmin in evaluating 
the selection function. We adopt that result here: we include in the analysis all W91PP objects down 
to the UGC limit, but set the formal diameter limit for evaluating the selection function to 1.15 arcmin. 
As with the MAT sample, we require log(a/6) > 0.1. The total number of W91PP galaxies thus included 
is 247. 

The selection criteria for the CF sample also included a photographic magnitude limit of 15.5. The 
expressions derived by Willick (1994) for dealing with this two-limit case exactly are unfortunately not 
analytic, and therefore are unsuitable for VELMOD. We thus decided to cut the CF sample at a larger 
diameter so that the magnitude limit would be relatively unimportant, and then use the one-catalog 
selection function corresponding to this larger diameter. Specifically, we include only those CF objects 
with UGC diameters > 1.5 arcmin in the VELMOD sample, and use a value of 1.6 arcmin in computing 
the CF selection function to account for residual incompleteness near the limit. As before, we also 
require log(a/6) > 0.1, and the total number of CF galaxies included in the VELMOD analysis is 170. 

Our method of assigning diameter limits in the W91PP, CF, and MAT for computing sample selection 
functions is far from satisfactory^. At distances ^ 5000 km s" 1 selection bias becomes an important effect 
for the forward relation. The treatment of selection effects in VELMOD is in principle rigorous, but 



is correct only to the degree that sample selection is properly characterized. We argued in § 2.2 that 



we can test our sensitivity to these problems by carrying out our analysis using both the forward and 



3 The majority of galaxies in A82 are close enough that the selection biases are not a serious issue. 
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inverse relations. As we shall see, the results for the two are in excellent agreement, which implies that 
our modeling of selection is adequate. 

The total number of galaxies that enter into the current analysis is 1876. The TF subsamples, their 
selection criteria, and the number of objects involved in each are summarized in Table 1. As discussed 
in Paper I, the cluster samples in the Mark III Catalog, HMCL and W91CL (cf. Willick et al. 1997a), 
are not suitable for the VELMOD approach, which is tailored to field galaxies, and we do not include 
those samples here. We also have elected not to include the elliptical galaxy portion of the Mark III 
catalog. 



4. Treatment of the Quadrupole 

In Paper I we presented evidence of systematic residuals from the IRAS-predicted velocity field, 
which could be modeled as a velocity quadrupole of the form uq = Vqt, where Vq was a traceless, 
symmetric 3x3 matrix. We argued that the probable cause of this quadrupole residual was differences 
between the true and measured density fields due to shot noise and the smoothing process, at distances 
> 3000 kms" 1 . 

We initially treated the five independent components of Vq as free parameters at each value of (3j in 
our Paper I analysis. In our final VELMOD run, however, we held the quadrupole fixed at the values 
of the five components obtained by averaging their maximum likelihood values at each (3j, so that our 
quadrupole residual would not fit out the /^-dependent part of the quadrupole already present in the 
IRAS-predicted velocity field. The final values of these quadrupole components are given in Paper I, 
Table 2, and a map of the overall quadrupole flow is presented in Figure 4 of Paper I. Its rms amplitude 
averaged over the sky is 3.3% of Hubble flow. 

The Paper I quadrupole increases linearly with distance, which is the expected signature of a 
quadrupole generated at distances beyond the sampled region. However, there is no reason to be- 
lieve that the linear quadrupole extends beyond 3000 kms -1 . A substantial fraction of the Paper I 
quadrupole is generated by mass density determination errors at distances 3000 ;$ r <, 12000 kins -1 (Pa- 
per I, Appendix B). Such errors will give rise to a linear quadrupole only at r^,3000 kms" 1 . In the 
region coincident with the mass determination errors, the velocity residual will not have a quadrupole 
form at all, and in fact will not be expressible as a divergence-free flow. Only at distances beyond the 
region of dominant mass determination errors will a divergence-free velocity residual reappear, now with 
an r~ 4 dependence rather than a linear one (cf. Jackson 1976, equation 3.70). 

For this paper we adopt the simplest model consistent with both our Paper I result and the above 
considerations. We assume that the radial component of the residual velocity field is given by 

Mo(r) = — — r-r, (6) 

where r = r/r. In equation @, Vq is the same traceless, symmetric 3x3 matrix as was derived in 
Paper I. However, we have introduced a new quantity, Rq, that parameterizes the cutoff scale of the 
linear quadrupole. For r <C Rq, we recover the Paper I quadrupole exactly. For r 3> Rq, we obtain 
the r -4 quadrupole expected at large distances. The transition between them is smooth but rapid, so 
there is only a small region, with r Rq, in which uq(t) does not behave like a quadrupole. This 
is a desirable feature, for it minimizes the volume in which our residual velocity field has divergence. 
We will determine the value of Rq through maximum likelihood in § 5.3, and justify our modeling of 
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the quadrupole a posteriori in § |7], where we show that the IRAS velocity field, with the quadrupole 
included, gives an acceptable fit to the TF data. 



5. Results 



In this section we present the main results of applying VELMOD to the 1876-galaxy, czlg < 
7500 kins -1 subsample described in § [| In § |5.1| , we search for a luminosity dependence to the TF 
scatter. We show the results for (3j without allowing for an external quadrupole in § |5.2| . In § |5.3| we 
find the value of Rq to use in the quadrupole formula, equation ||. We then use this quadrupole in 
an analysis of the small-scale velocity dispersion in § where we give our final results for f3j. The 
robustness of this result to subsample is discussed in § |5.5|, and to smoothing scale in § ^6. 



5.1. Establishing the luminosity/velocity-width dependence of the TF scatter 

Giovanelli et al. (1997) and Willick et al. (1997a) have pointed out that in some samples, the TF 
scatter is a function of luminosity. In Paper I, we modeled the velocity-width dependence of the forward 
TF scatter as having the form <7tf(?7) = °o — 9fVi an d adopted the Willick et al. (1997a) values of 
gf = 0.14 for A82 and gf = 0.33 for MAT. Here we similarly model the luminosity-dependence of the 
inverse TF scatter by a v (M) = a V fi + 9i(M — M), where M is the mean absolute magnitude for the 
sample^], but determined the gf and gi through maximum likelihood, as follows. 

First, we ran a preliminary set of VELMOD runs, both forward and inverse, in which the gf and the 
gi were treated as free parameters at each value of (3j. These runs demonstrated that there was negligible 
cross-talk between the gf (g^) and any other parameter of interest, in particular, We thus used these 
preliminary runs to establish their values, and then held them fixed for all subsequent VELMOD runs. 
The preliminary runs employed the simplest velocity models: no quadrupole, wlg fixed at its Paper I, 
no quadrupole value, and a v fixed at 150 kins" 1 without allowance for a density dependence (cf. § |5^ ). 

We imposed an additional constraint on the gf and g^. The TF relation implies that 

_ da^ _ d(a TF /b) 1 da TF _ gj_ 
9i ~ dM~ d(brj) b 2 dr] ~ 6 2 ' U 

where b is the forward TF slope for the sample in question. We can regard gf and gi as well-determined 
from the data to the degree that this relation holds. We found that the A82 and MAT samples satisfied 
equation (|7j), and thus adopted the values of gf and gi determined from the preliminary VELMOD 
runs for those samples. For W91PP, however, gf and gi obtained from those runs had opposite signs 
and were thus inconsistent with equation (]?]). We interpret this to mean that there is no significant 
luminosity or velocity width dependence of the W91PP TF scatter, a conclusion also reached by Willick 
et al. (1997a). We thus set g f = g t = for W91PP. 

For CF, the preliminary gf and gi were both positive, but was significantly smaller than b~ 2 gf. CF 
uses optical widths, as does MAT. MAT shows the strongest signal of a luminosity-dependent scatter, 
which we conjecture is a consequence of optically-measured widths. We thus assigned CF the same value 
of gf as was found by maximum likelihood for MAT. For the CF gi, we took the mean of its maximum 
likelihood value and the value inferred from equation (|7|) given the adopted value of gf. 



6 The absolute magnitudes were calculated for this purpose as M = m— 51ogczLG, so that they are independent of /3/. 
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TABLE 2 

Luminosity /Velocity-width 
Dependence of TF Scatter 



Sample 


9f 


9i 


A82 


-0.24 


-0.0021 


MAT 


0.35 


0.0055 


W91PP a 


0.00 


0.0000 


CF 


0.35 


0.0030 



Table 2: Notes: (a) No significant luminosity/ width dependence of scatter was detected. 

We summarize the results of this exercise in Table |2|. Column (1) gives the sample name, while 
columns (2) and (3) list the adopted values of gf and gi respectively. Note that the MAT value of gj 
is very close to that found by Willick et al. (1997a). This is not surprising, as MAT is the only sample 
for which the luminosity dependence of scatter has a strong, unambiguous signal. The A82 value of 
gf, however, not only differs from the Willick et al. (1997a) value, but is negative, signifying a scatter 
increase with increasing luminosity. The physical reason for this is unclear, but the consistency between 
gf and gi for A82 suggests that the effect is real. However, these choices have no meaningful effect on 
the derived values of fli (as mentioned above) or any other important quantity discussed in this paper. 
Values of o"tf or a n quoted later in this paper refer to their values at r\ = and M = M respectively. 

5.2. No-quadrupole results 

After adopting the values of g/ and gi in Table 2, we reran forward and inverse VELMOD with no 
quadrupole and a v fixed at 150 kms" 1 . For these runs, we also fixed the LG random velocity vector 
w LG at the value determined in Paper I for the no-quadrupole case. 

The results are shown in Figure |T] for the forward (left panel) and inverse (right panel) TF relations^. 
The full likelihood versus f3i curves are quite similar for the forward and inverse TF relations. In 
particular, the maximum likelihood values of /?/ differ by only 0.01, which is insignificant given that the 
1 a error in (3i is 0.04. As we shall see, forward and inverse give essentially identical results for /3j for 
all VELMOD runs. Agreement between the forward and inverse results means that our approximate 
treatment of the selection and luminosity functions have no meaningful effect on f3i (see the discussion 
in§P). 

Finally, note that the value of f3i obtained here for the no-quadrupole case is very close to the value 
of Pi = 0.56 obtained in Paper I for the no-quadrupole case. Thus, more than doubling the number of 
sample objects and extending the redshift limit from 3000 to 7500 kms" 1 has had essentially no effect 
on ftj (other than shrinking the error bar) . We will see below that the same is true when the quadrupole 
is included. 



7 Thc absolute values of the forward and inverse likelihood statistics are quite different because the former derives from 
a probability density in apparent magnitude, the latter from a probability density in the width parameter rj. 
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Fig. 1. — The VELMOD likelihood statistic £ as a function of Pi for the no-quadrupole case. Results are 
shown for both the forward (left panel) and inverse (right panel) TF relations. In both cases the velocity noise 
a v has been held fixed at 150 kms -1 and the LG random velocity vector wlg has been fixed at its Paper I, 
no-quadrupole value. The dotted curves are the cubic fits to the C(fti) points. The maximum likelihood values 
of j3i are indicated on the plots as (3 m \. 



5.3. Determining the quadrupole cutoff scale 

As described in § U we adopted a quadrupole velocity residual, equation (|6|), that agrees with the Paper 
I quadrupole at small distances, but changes smoothly from a linear to an r~ 4 quadrupole for r > Rq. 
To determine the value of Rq we carried out a series of VELMOD runs, both forward and inverse, with 
values of Rq ranging from 100 kms -1 (which essentially means no quadrupole) to 15, 000 kms -1 (which 
amounts to the Paper 1 quadrupole throughout the sample volume). In each of these runs, Rq was held 
fixed, as were a v at 150 km s" 1 and wlg at its best-fit Paper I value when the quadrupole was included. 
Only j3i and the 12 TF parameters were varied in each run of a given Rq. 

Figures and |3] show the results of these runs for the forward and inverse relations respectively. In 
each figure, the upper panel shows £ versus Rq, while the lower panel shows the maximum likelihood 
value of Pi versus Rq. We show the results only out to Rq = 10,000 kms -1 , as for larger Rq neither 
Pi nor C changed appreciably. 

Note that Pi is very insensitive to Rq] over the entire range of Rq considered, Pi changes only by ~ 
0.03, or less than 1 a. There is a well-defined likelihood maximum (minimum of C) at Rq = 3000 km s _1 
for the inverse case and at Rq = 4000 kms -1 for the forward case. Note that Rq < 2000 kms -1 and 
Rq > 5000 kms -1 are strongly disfavored in both cases, while Rq = 3500 kms -1 is consistent with 
both, so we adopt the latter value for the remainder of the paper. The maximum likelihood values of 
Pi are very close to 0.50 for this value of Rq, for both forward and inverse VELMOD. 

The value of Rq = 3500 kms -1 signifies that the Paper I quadrupole cuts off strongly beyond this 
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Fig. 2. — The best forward TF likelihood £f or w (upper panel), and the corresponding maximum likelihood 
value of (3i (lower panel) plotted as a function of the quadrupole scale Rq (cf. equation Q). For these runs 
the velocity noise was held fixed at a v = 150 kms -1 and the LG random velocity vector wlg at its Paper I, 
quadrupole value. Note the strong likelihood maximum for Rq = 3000-4000 kms -1 . The corresponding value 
of Pi is 0.49. 

distance. We will see in § that the resulting velocity field is an adequate fit to the data. We can therefore 
conclude that the IRAS versus true mass differences arising from the smoothing/filtering procedure that 
dominate the velocity prediction errors are concentrated in the range ~ 3000-5000 kms -1 . 

At Rq = 100 kms" 1 , corresponding to essentially no quadrupole, we find fti = 0.48 for the forward 
relation, which differs from the no-quadrupole value of 0.53 found in § |5.2j . These differ because different 
values of wlg were used: In § |5]2| we fixed wlg at the Paper I, no-quadrupole best-fit value, while here, 
we used the Paper I quadrupole best-fit value. The values of £ in the two cases are similar, however, 
implying that we have limited sensitivity to wlg in the likelihood analysis. There is moderate covariance 
between wlg and j3i when a quadrupole is not used to describe the local flow field. With the quadrupole 
added, however, this covariance is much reduced. Stated another way, when the quadrupole is modeled, 
wlg is both smaller in amplitude and better determined. 
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Fig. 3. — As in the previous figure, for the inverse TF relation. 

5.4. The small-scale velocity dispersion 

In the VELMOD runs described up to now, the small-scale velocity dispersion a v has been fixed at 
150 kms -1 , a useful round number with which to establish the values of gf, g^, and Rq. Having done 
so, we ran a series of VELMOD runs for a range of fixed values of a v , with Rq fixed at 3500 kins" 1 and 
Wlg fixed at its Paper I, quadrupole value. In each run, j3i and the 12 TF parameters were varied to 
maximize likelihood. Figures § and |^ show the results for the forward and inverse relations respectively. 
In each case, we plot the maximum likelihood values of f3i and the corresponding C(/3i) versus a v . 

There is a weak systematic variation of j3i with a v , amounting to less than 0.05 over the full range 
of a v considered. The likelihood reaches a clear maximum at a v = 150 ± 20 kms" 1 for the forward 
relation, and 130 ± 20 kms -1 for the inverse relation. Over the la errorbar, the maximum likelihood 
values of /3j only vary by 0.01, much less than the statistical error on this quantity. Thus there is little 
covariance between f3j and o v . 

The value of a v found here is consistent with the maximum likelihood value a v = 125 ± 20 kms -1 
found in Paper I. This is reassuring, though not surprising; as discussed in Paper I, a v is primarily 
determined at small distances, < 3000 kms -1 , where its effect on the overall variance is comparable to 
that of the TF scatter. 
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Fig. 4. — The best forward TF likelihood £f or w (upper panel), and the corresponding maximum likelihood 
value of 0i (lower panel) plotted as a function of the small-scale velocity dispersion a v . For these runs the 
quadrupole scale was fixed at Rq = 3500 kms , and the LG random velocity vector at its Paper I, quadrupole 
value. The likelihood statistic yields a v = 150 ± 20 kms -1 . Within this favored range 0i varies trivially, from 
0.48-0.49. 

5. 4-1- Density-dependence of a v 

Strauss, Ostriker, & Cen (1998) and Kepner, Summers, & Strauss (1997) showed that the small- 
scale velocity dispersion is an increasing function of local density. In Paper I, we chose to neglect such 
variation, the only exception being our "collapsing" of 20 Virgo cluster galaxies by assigning them 
redshifts equal to the cluster mean (cf. Paper I, §4.3). For this paper we attempt to detect a density- 
dependence of a v through the likelihood analysis. We adopt a model of the form 

<7v(8g) = <7v,l +fs{Sg ~ 1) , (8) 

where 5 g is at the same smoothing as was assumed for the IRAS velocity field calculation. We take 
5 g = 1, rather than 5 g = 0, as the zero point for our model because most TF sample objects lie in 
relatively high-density environments (the mean value of S g for the full TF sample is ~ 0.8). 

We carried out a series of forward and inverse VELMOD runs for a range of values of o~ v> i. In each 
run Pi, the twelve TF parameters, and f$ were treated as free parameters. We continued to hold Rq 
fixed at 3500 kms -1 . The results of this exercise are shown in Figure |6] for the forward relation, which 
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Fig. 5. — As in the previous figure, for the inverse TF relation. 



plots the maximum likelihood values of fii, and the corresponding values of fg and Cf OTW as a function 
of a Vi \. Allowing for a density-dependent velocity dispersion has no significant effect on our derived 
which remains very close to 0.50 near the minimum of Cf OTW . 

The best likelihood is achieved for a Vj i ~ 140 kms -1 , very similar to the value of the invariant a v for 
which likelihood was maximized (compare with Figure ^). For a V) i in the range favored by the likelihood 
statistic, 140 ± 25 kms -1 , fs is remarkably constant at 33-34 kms -1 . The minimum value of Cf ovw in 
Figure ^ is 5.5 points smaller than its minimum value for an invariant a v , corresponding to an increased 
likelihood of the fit by a factor of e*- 5 ' 5-1 ^ 2 ~ 9.5, a 2.1<r result. For the inverse relation (not shown) 
fs = 36 kms -1 when the best likelihood is achieved for a V) i = 120 ± 25 kms -1 , and the best likelihood 
is greater than for the invariant a v case by a factor of ~25, a 2.5 a result. We have thus detected a 
significant variation of velocity dispersion with density. To a good approximation we may summarize 
these results (now normalizing to 5 g = 0) as a v = [100 ±25 + 35 S g ] kms -1 . 

The value of a v for galaxies in a mean density environment is very small, consistent with the conclu- 
sions of Paper I, Davis, Miller, & White (1997), Strauss et a/. (1998), and papers referenced therein. 
The quantity a v is the quadrature sum of IRAS error and true velocity noise (cf. Paper I, § 3.2). We 
estimated the former to be ~ 70 kms -1 in Paper I, so the true 1-D velocity noise is only about 50-70 
kms -1 in mean-density environments. The flow field of galaxies is remarkably cold. 

In Figure [7| we plot the forward and inverse likelihood statistics versus f3j. The plots are done for 
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Fig. 6. — The results of a series of forward VELMOD runs in which the small-scale velocity dispersion is treated 
as a linear function of galaxy density (cf. equation Q). The bottom panel shows the maximum likelihood value 
of Pi, the middle panel the corresponding value of fs, and the upper panel the corresponding value of the 
likelihood statistic as a function of cr„,i- The Paper I quadrupole with Rq = 3500 kms" 1 was used. 



a Vy i = 130 kms" 1 in the forward case, and a Vt i = 110 kms -1 for the inverse case. As in Figure [l], the 
forward and inverse curves are almost identical, and the resultant maximum likelihood values of fli are 
the same to within 0.01 (see Table |3|). This tells us that any errors we may have made in modeling 
sample selection and luminosity functions have had little or no effect on the quantities of interest. 



5.5. Breakdown by sample and redshift 

We test the robustness of our results by computing the maximum likelihood f3j for different TF 
subsamples and redshift ranges. The is done in Table for our favored forward and inverse runs. 
Each of the four TF subsamples, for both the forward and inverse TF relations, produces a maximum 
likelihood f3i consistent with one another and with the global value of 0.50. Similarly, the maximum 
likelihood j3i for objects in each of five redshift bins are statistically consistent with one another. The 
last redshift bin gives a value of f3j somewhat higher than the others, but the error bar is larger, and it 
is still consistent. Thus, there is no significant trend with redshift. This consistency among sample and 
redshift range enhances our confidence in our global value of j3i. Note that lower-redshift galaxies give 
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Fig. 7. — The VELMOD likelihood statistic £ as a function of (3i for the forward TF (left panel) and inverse TF 
(right panel) relations. The maximum likelihood values of /3/ are indicated on the plots. For these VELMOD 
runs, the 300 kms _1 -smoothed IRAS velocity field, the Paper I quadrupole with Rq = 3500 kms" 1 , and the LG 
random velocity vector wlg fixed at its Paper I value were used. A density-dependent velocity dispersion (cf. 
equation ||) was used in these runs, with f$ treated as a free parameter at each /?/. A value of f$ « 35 kms -1 
was obtained for fii = 0.5 in both cases. For the forward TF run the results for a v \ = 130 kms -1 are plotted, 
while for the inverse TF run the results for a V) i = 110 kms -1 are shown. 

more leverage per object on f3j than do higher redshift galaxies. The reasons for this were discussed in 
Paper I, §4.5. The W91PP sample yields the weakest constraints on because of the relatively small 
volume it probes. Although there are fewer CF than W91PP galaxies, they yield a stronger constraint 
on f3i because the CF sample has wider sky coverage. 

5.6. Results for 500 kms 1 smoothing 

Our Paper I tests with mock TF and IRAS catalogs showed that VELMOD returned unbiased esti- 
mates of fii when a 300 kms -1 Gaussian smoothing scale was used in the IRAS velocity predictions. 
We also tested a 500 kms -1 smoothing scale and found that it produced estimates of /3j biased ~25% 
high. For the real data, 500 kms" 1 smoothing produced a maximum likelihood f3j about 15% higher 
than the 300 kms -1 value (cf. Paper I, § 4.6). 

The results of applying VELMOD to the expanded sample using 500 kms _1 -smoothed IRAS velocity 
predictions are shown in Figure in which likelihood for the forward and inverse TF relations is 
plotted versus These runs are carried out using the values of the quadrupole parameters and wlg 
obtained from the Paper I 500 kms" 1 run, but again using the modified quadrupole of equation (H) 
with Rq = 3500 kms" 1 . We again allow for a density-dependent value of <r v ; we show the results for 
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TABLE 3: Breakdown by Sample and Redshift 



Subsample 


f3i (forward) 


(3i (inverse) 


N 


A82 


0.477 ±0.062 


0.486 ±0.061 


300 


MAT 


0.518 ±0.052 


0.533 ±0.052 


1159 


W91PP 


0.411 ±0.159 


0.386 ±0.148 


247 


CF 


0.488 ±0.107 


0.478 ±0.117 


170 


cz < 1500 kms -1 


0.515 ±0.059 


0.510 ±0.056 


327 


1500 < cz < 3000 kms -1 


0.542 ±0.065 


0.532 ±0.066 


564 


3000 < cz < 4500 kms -1 


0.428 ±0.084 


0.473 ±0.088 


370 


4500 < cz < 6000 kms -1 


0.376 ±0.096 


0.381 ±0.094 


422 


6000 < cz < 7500 kms -1 


0.594 ±0.173 


0.734 ±0.197 


193 


Overall 


0.495 ±0.037 


0.503 ±0.036 


1876 



Table 3: Notes: Results are given for the preferred VELMOD runs: 300 kms _1 -smoothed IRAS 
predicted velocity field; density-dependent velocity dispersion with a v> i = 130 kms -1 (forward) and 
a v> i = 110 kms -1 (inverse); Paper I quadrupole with Rq = 3500 kms -1 , and corresponding value of 
wlg- The subsample /Vs were calculated using the TF and velocity parameters obtained from the 
full-sample run; these parameters were not solved for separately for each subsample. 
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Fig. 8. — As in the previous figure, now using the 500 kms -1 -smoothed IRAS velocity field predictions. The 
maximum likelihood values of (3i differ by very little from those of the previous figure. 



<7 V) i = 150 kms -1 , which maximizes likelihood at 500 kms -1 smoothing. 

The 500 kms -1 maximum likelihood estimates of /3i differ little from those obtained at 300 kms -1 
smoothing. Averaging the forward and inverse results, we find /3i = 0.52 ± 0.05, only 4% higher than 
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TABLE 4: VELMOD and Mark III TF Relations' 





forward 






inverse 




Sample 


A 


b 


0"TF 


D 


e 


a v 


A82 (VELMOD) 


-5.96 


10.44 


0.45 


-5.96 


0.0879 


0.042 


A82 (Mark III) 


-5.94 


10.29 


0.47 


-5.98 


0.0893 


0.043 


MAT (VELMOD) 


-5.80 


7.16 


0.43 


-6.00 


0.1282 


0.063 


MAT (Mark III) 


-5.79 


6.80 


0.43 


-5.96 


0.1328 


0.059 


W91PP (VELMOD) 


-4.09 


7.14 


0.40 


-4.13 


0.1217 


0.052 


W91PP (Mark III) 


-4.28 


7.12 


0.38 


-4.32 


0.1244 


0.049 


CF (VELMOD) 


-4.00 


8.41 


0.48 


-3.97 


0.0948 


0.049 


CF (Mark III) 


-4.22 


7.73 


0.38 


-4.27 


0.1190 


0.047 



Table 4: Notes: (a) Comparison of the VELMOD and Mark III TF calibrations for the four subsamples 
used in the VELMOD analysis. The typical 1 a errors for both calibrations are: 5 A = SD ~ 0.03; 
S b /b = 5e/e = 0.03; 5a TF = e~ l 5a v = 0.02. 

our 300 kms -1 result. In the VELMOD analysis, the TF data are not smoothed, and therefore we chose 
a small smoothing scale for the IRAS density field in order to model the velocity field in as much detail 
as possible. The fact that Pi, and, more significantly, the best values of £{ orw and C imn are essentially 
unchanged when the smoothing is increased to 500 kms -1 , says that density fluctuations on scales 
between 300 and 500 kms -1 contribute little to the velocity field. That is, there is little small-scale 
power both in the true velocity field, and in the IRAS-predicted velocity field (i.e., the gravity field). 
The simulations used in the mock catalogs in Paper I do have a substantial amount of small-scale power 
in the density field, and this is presumably the reason that they yielded a substantially biased estimate 
of j3i, and a far worse value of C, with 500 kms -1 smoothing. Because of the lack of small-scale power 
in the velocity field, the agreement between our 300 and 500 kms -1 results for /?/ does not shed light 
on the question of whether biasing is scale-dependent. 

Our 500 kms -1 runs also detect an increase in the small-scale velocity dispersion with density. For 
the forward run we find fs = 30 kms -1 , similar to the 300 kms -1 value. However, for the inverse run we 
find fs = 60 kms -1 . We would expect fs to rise with smoothing scale because the density contrasts are 
generally smaller with larger smoothing. The inverse result confirms this expectation but the forward 
does not; we do not understand the reason for this difference. 

6. VELMOD versus Mark III Catalog TF Calibrations 

An important feature of VELMOD is that the TF relations for the various samples are determined 
by maximizing likelihood at each /?/. The correct TF relations are those obtained for the maximum 
likelihood value of j3j, which we have found to be 0.50 with small uncertainty. In Table [| we give the 
parameters of the forward and inverse TF relations obtained from our favored VELMOD runs, i.e. 300 
kms -1 smoothing, quadrupole with Rq = 3500 kms -1 , and density- dependent velocity dispersion with 
a Vt i = 130 kms -1 (forward) and a Vy \ = 110 kms -1 (inverse). Columns (1), (2), and (3) list the forward 
TF parameters A, b, and o"tf, while columns (4), (5), and (6) list the inverse TF parameters D, e, 
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Fig. 9. — Differences in the radial peculiar velocity inferred from the Mark III and VELMOD forward TF 
calibrations, plotted as a function of Local Group redshift for the A82 (left panel) and MAT (right panel) 
samples. 

and cr^f]. Also given in Table 3 are the values of these parameters that went into the Mark III Catalog 
(Willick et al. 1997a). 

The slopes and scatters of the VELMOD and Mark III TF relations are in good agreement overall. 
The VELMOD MAT TF slope is higher, but by less than 2 a, as we discussed in Paper I, §4.7. The CF 
slope is higher than its Mark III value, as is its scatter. This is not surprising, because in this paper we 
have treated CF as a fully independent sample, whereas the Mark III calibration procedure (Willick et 
al. 1996) assumed that CF had the same TF relation as the Willick (1991) cluster sample, W91CL, up 
to a slight zero-point adjustment. 

More important, there are substantial zero-point differences between the VELMOD and Mark III 
calibrations. While the VELMOD and Mark III TF zero points of A82 and MAT are in good agreement, 
those of W91PP and CF differ by about 0.2 mag, for both the forward and inverse TF relations. This 
difference is much greater than the expected errors of ~0.03 mag in either procedure. Because the 
difference manifests itself for only two of four samples, it cannot arise from a global zero point error in 
either the Mark III or the VELMOD calibration procedure. 

Figures |9| and |It] show how these differences in the TF parameters translate into peculiar velocity 
differences. The differences between the Mark III and VELMOD peculiar velocities inferred from the 
forward TF relation are plotted as a function of LG redshift for each of the four samples. The plots 
would appear substantially the same if we used inverse TF distances. We do not apply Malmquist 



The scatters are given for r\ — and for the mean absolute magnitude in each sample; cf. § 5.1 
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Fig. 10. — As in the previous figure, for the W91PP (left panel) and CF (right panel) samples. 



bias corrections, which would accentuate the differences between the VELMOD and Mark III velocities. 
Thus, the TF scatters have no effect on the diagrams. 

For A82 there is no meaningful difference between the Mark III and VELMOD inferred peculiar 
velocities. For MAT, there is a slight trend, but the mean differences are everywhere less than ~ 100 
kms -1 , except at the outer edge (cz ~ 6000 kms -1 ) of the sample. However, for W91PP and CF the 
differences are substantial. In each case, the Mark III velocities are more negative by 200-400 kms -1 . 
In the case of W91PP, the differences are even larger beyond 6000 kms -1 .[] 

This systematic difference between the Mark III and VELMOD TF calibrations has a strong effect 
on the inferred bulk flow from the Mark III data (cf. Courteau et al. 1993; Dekel 1994; Postman 1995; 
Strauss 1997, for discussions). The W91PP and CF samples dominate the Northern sky away from 
the Local Supercluster. W91PP in particular samples the Perseus-Pisces (PP) supercluster, centered at 
/ rs 120° , b w —30° . As measured by the Mark III TF calibrations, the PP region is seen as having large, 
negative radial peculiar velocities in the microwave background frame (e.g., Courteau et al. 1993). This, 
along with outflowing velocities in the Great Attractor region (traced mainly by the MAT sample), is 
why measurements of the bulk flow within 6000 kms -1 from the Mark III data have yielded values in 
the range 400-500 kms -1 . However, IRAS does not predict strong infall of the PP supercluster region, 
unless (3j is ;S0.2. Since the VELMOD TF calibrations reflect the IRAS velocity field, they adjust to 
produce little infall of PP, and thus a much smaller bulk flow, than do the Mark III calibrations. 



9 For W91PP the trend is essentially linear with redshift, and has small scatter, whereas for CF, there is larger scatter 
and the velocity difference levels off at large redshift. This is because for W91PP the calibration difference involves only 
the TF zero point, while for CF both zero point and slope differences are present. The TF slope difference also explains 
why the MAT diagram exhibits a much larger scatter than the A82 diagram. 
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Another way to state the problem is as follows. The Mark III TF zero points were set by asking for 
agreement in distances for galaxies in overlapping datasets; the full-sky cluster sample of Han & Mould 
(1992; HMCL) was the backbone that tied the sky together (cf. Willick et al. 1995, 1996, 1997a). If these 
calibrations are indeed correct, then the VELMOD calibrations are not, and it follows that the IRAS 
redshift survey does not correctly predict the peculiar velocity field. In fact, this was the conclusion 
reached by Davis et al. (1996), whose ITF analysis made use of the Mark III zeropointing procedure 
even though it did not use the Mark III distances directly. If the IRAS velocity field predictions are 
correct, as we have assumed in this paper, then so are the VELMOD TF calibrations and our maximum 
likelihood estimate of f3[. However, in that case the Mark III TF calibrations are incorrect, and the Mark 
III Catalog contains erroneous distances for the W91PP and CF samples — and by extension, for the 
HMCL, W91CL, and elliptical galaxy samples as well. It would then follow that the POTENT peculiar 
velocity and density fields, which are based on the Mark III distances and were used in the POTIRAS 
determination of f3j = 0.89 ± 0.12, contain systematic errors. A self-consistent picture would require 
that the VELMOD TF calibrations, required by the IRAS velocity fields, also be used to produce the 
POTENT velocity and density maps to estimate This has not yet been done. 

One can ask whether the VELMOD TF calibrations agree better with the Mark III calibrations for 
some value of f3j other than 0.5. In fact, for (3j = 0.1, the VELMOD W91PP zero point agrees with 
that of Mark III, for both the forward and inverse TF relations. However, for /3j = 0.1 the VELMOD 
TF zero point for CF is even farther from its Mark III value than it is for fii = 0.5. For (3i ~ 1, the 
CF zero point is closer to its Mark III value, but the W91PP zero point diverges drastically from Mark 
III. Also, for very low or very high (3i we lose the good agreement between the VELMOD and Mark 
III A82 and MAT TF zero points. Thus there is no value of j3i at which the VELMOD and Mark III 
calibrations are in overall agreement. 

The question of which set of TF calibrations is correct must ultimately be decided by improved TF 
data. The problem has arisen because there is no reliable way to tie together the disjoint Southern 
(MAT) and Northern (CF and W91PP) sky TF data sets that constitute the Mark III field spirals. A82 
spans the two hemispheres but is dominated by nearby galaxies and has little overlap with the Northern 
sky samples. The HMCL sample was thought to provide the needed overlap, but its uniformity across the 
sky has been called into question by the calibration disrepancies. What is needed are homogeneous TF 
data that cover the celestial sphere. In collaboration with S. Courteau, M. Postman, and D. Schlegel, we 
have obtained uniform TF data for ~ 300 galaxies isotropically distributed in the spherical shell defined 
by 4500^,02^7000 kms -1 . Reduction of these data are under way, and results are expected by late 
1998. Comparison of these uniform TF data with the Mark III data will allow a definitive resolution of 
the calibration problem. 

Finally, we note that adopting the Mark III TF calibrations has relatively little effect on the maximum 
likelihood /?/ obtained from VELMOD. With the 300 km s~ ^smoothed IRAS plus quadrupole velocity 
model, we obtain f3i = 0.44 (forward) and (3i = 0.45 (inverse) when the TF parameters for all four 
samples are fixed to their Mark III values as given in Table |4j. For the no-quadrupole model we obtain 
/3j = 0.50 (forward) and f3j = 0.51 (inverse). The likelihoods obtained from these VELMOD runs are, of 
course, much worse (by ~ 100 units in C) than for our preferred runs in which the twelve TF parameters 
are free. Thus, while the TF calibration problem is crucial for the match of the IRAS velocity field to 
the TF data, as we discuss in the next section, it is secondary for the determination of 
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7. The Goodness of Fit of the IRAS Velocity Field 

Although VELMOD does not produce a picture of the TF velocity field, we can nonetheless use it to 
visualize how well the TF data fit the IRAS velocity predictions. We do so by converting the VELMOD 
apparent magnitude m (forward) or velocity width parameter rj (inverse) residuals into smoothed radial 
peculiar velocity residuals with respect to IRAS, as described in Paper I, § 5.1. The VELMOD residuals 
also enable us to measure the goodness of fit of the velocity model, as we describe below. The smoothed 
peculiar velocity residual is given by equation 24 of Paper I: 



5u s , = dj 



l-/aO a2(<5 ^ xAmi) l , (9) 



which we repeat here because of a typographical error in Paper I; see Paper I for the definition of the 
various symbols in this equation. 

Figures [TT|, [T2|, and [13] show sky maps of these velocity residuals for f3i = 0.5, (3j = 0.1, and (3j = 1.0 



respectively. In each case, the results are based on forward TF residuals from our preferred 300 kms -1 
smoothing run (see the notes to Table 3). Open symbols represent negative velocity residuals (i.e., the TF 
distance to the object is greater than that predicted by IRAS); starred symbols represent positive velocity 
residuals. The Gaussian smoothing scale for the maps is given by 250 [1 + (cz LG /2500) 2 ] 1/2 kms" 1 . Thus, 
the smoothing radius varies from 250 kms -1 nearby to ~ 750 kms -1 at the edge of the sample. This 
smoothing imposes a coherence scale of ~ 15-25° on the results; patches this size with similar velocity 
residuals are to be expected in the maps from the smoothing alone, while any coherence seen on much 
larger scales represents a real error in the model. Points are plotted only for galaxies which have enough 
near neighbors to allow an adequate smoothing; this is why there are few galaxies represented in the 
Northern Galactic Cap at cz > 2500 kms -1 , where the sampling is very dilute. Such points, if plotted, 
would exhibit large velocity residuals due solely to TF scatter and would not help us assess the quality 
of the fit. 

Inspection of these maps shows clearly why f3j = 0.5 is the best fit. Although there is some real excess 
coherence to the residuals (we discuss this further below), the coherent velocity levels are generally at a 
low level ( ;$ 250 kms -1 ). There are many alternating regions of positive and negative residuals, showing 
that globally at least the residual map is fairly incoherent. This is what is required of a good fit. There 
are no extended regions where the velocity residuals are consistently greater than 300 kms -1 . This is a 
qualitative indication that the IRAS plus quadrupole velocity field model fits the major features of the 
actual velocity field. As in Paper I, coherent residuals are present when the quadrupole is not modeled. 
That being said, with Rq = 3500 kms -1 , the quadrupole contribution at cz ^5000 kms -1 is negligible. 
Thus, the good agreement on very large scales is due to the IRAS velocity field alone, giving a posteriori 
confirmation of our quadrupole model. 

The residual maps produced at (3j = 0.1 and at /3j = 1.0, on the other hand, show considerable 
coherence. Moreover, the amplitude of the velocity residuals in these regions is often large, ^ 300 km s -1 . 
Low and high f3j are clearly worse fits to the TF data than is (3j = 0.5. The maps, then, confirm what 
the likelihood analysis is telling us. It is important to remember that the poor fit at low and high (3i is 
not a result of errors in the assumed TF relation, for the TF relations used were those preferred by the 
data at each The poor fit is a genuine reflection of the incorrectness of the IRAS velocity field for 
low and high 
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P , = 0.5 




Fig. 11. — Smoothed VELMOD velocity residuals plotted in Galactic coordinates, for f3i = 0.5. Open circles 
indicate objects inflowing relative to the velocity model, while starred symbols represent outflowing objects. 
The symbol size indicates the magnitude of the velocity residual, as coded at the lower right of each plot. 
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Fig. 12. — Same as the previous figure, but for f3j = 0.1. 
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Fig. 13. — Same as the previous figure, but for f3j = 1.0. 
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Fig. 14. — VELMOD residual autocorrelations, tp(r), plotted for /3j = 0.5, /3j = 0.1, and 0j = 1.0. Nonzero 
values of V ; (' r ) indicate coherence of the TF versus IRAS residuals on a spatial scale r. Such coherence is 
pronounced for /3j = 0.1 and Pi = 1.0, indicating a poor fit. It is insignificant for 0j = 0.5. 



We may quantify our visual impressions by means of the residual autocorrelation function ip(r), 



denned by equation (25) of Paper I. In Figure 14, we plot ip(r) for the three values of /3j represented 
in the previous figures. The plots show that for /?/ = 0.1 and /?/ = 1.0, significant excess correlation 
is evident on small and large scales. At /3j = 0.5, the ip(r) is consistent with zero on all scales. There 
is a small amount of positive correlation on scales ^2500 kms -1 for j3j = 0.5, consistent with the 
(low-amplitude) inflowing monopole residuals in the upper panel of Figure II. This may be indicative 
of a breakdown of the IRAS model at some level, but it is not highly significant, as we now show. 

A rigorous measure of the level of residual coherence comes through the use of the correlation \ 2 
statistic, x| ; defined by equation (26) of Paper I. We plot x\ versus Pi in Figure [15]. (The value for 
/3j = 0.1 is off-scale.) In Paper I we showed that this statistic had properties similar to that of a true 
X 2 statistic, but with a mean of 0.87 ± 0.06 per degree of freedom rather than unity. Its variance was 
consistent with that of a true \ 2 statistic. We indicate the expected value of x| (in this case, 63.5 for 
73 degrees of freedom) as a heavy solid line on the plot. The 1- and 3 o deviations from the expectation 
are indicated as dot-dashed and dashed lines, respectively. The quantity \\ reaches its minimum at the 
maximum likelihood value of The only other value of (5j for which x| is within 3a of the expectation 
value is (3i = 0.6. /?/ < 0.4 and (3i > 0.7 are ruled out at the > 3 a level. 
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Fig. 15. — The residual correlation statistic x|, plotted as a function of f3j. The heavy line shows the expectation 
value of this statistic; the dot-dashed and dashed lines show the 1 a and 3 a deviations respectively. 



8. Summary 

We have applied the VELMOD method to a TF sample drawn from the Mark III Catalog in order 
to estimate f3i = il os /bj, where bi is the linear biasing parameter for IRAS galaxies. The TF sample 
consists of 1876 galaxies, comprising nearly all Mark III field spirals to a limiting redshift of cz LG = 
7500 kms -1 . This analysis extends the one we presented in Paper I, which was limited to 838 galaxies 
with cz^g < 3000 kms -1 . As in Paper I, peculiar velocities were predicted from galaxy density contrasts 
obtained from the IRAS 1.2 Jy redshift survey (Fisher et al. 1995), under the assumption of linear 
gravitational instability theory and linear biasing. We developed an analytic approximation to the 
single-object VELMOD likelihoods, applicable to £ 75% of sample objects, which makes the code run 
3-4 times faster. 

We carried out the VELMOD analysis using both the forward and inverse forms of the TF relation. 
Consistency between the two is required to ensure that selection biases are unimportant. We found 
that the maximum likelihood values of as well as other important velocity parameters, were indeed 
statistically the same for both forms of the TF relation. In addition, we allowed the quadrupole velocity 
residual detected in Paper I to cut off smoothly beyond a radius, Rq, whose value we determined 
through likelihood maximization to be 3500 ± 1000 kms^ 1 . There is little covariance between Rq and 
Pi. We believe the quadrupole is real and readily accounted for (cf. Paper I, Appendix B). We may 
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summarize our results as /3j = 0.50 ± 0.04 ± 0.04, where the first errorbar is statistical and the second is 
systematic. This value is quoted for our favored model in which the IRAS densities are smoothed with 
a 300 kms -1 Gaussian, the small-scale velocity dispersion varies with density (see below), and in which 
the quadrupole, with Rq = 3500 kms -1 , is added to the IRAS-predicted velocity field. The systematic 
error is due to the quadrupole; if it is not valid to add it, we obtain f3i = 0.53 ± 0.04 (forward), or 
fii = 0.54 ± 0.04 (inverse). We also found that changing the IRAS smoothing scale from 300 to 500 
kms -1 does not significantly affect the derived value of f3j. This implies that there is little contribution 
to the velocity and gravity fields from fluctuations on scales between 300 and 500 kms" 1 . Further work 
needs to be done to quantify this, and to understand what effect our Wiener filter, which suppresses 
fluctuations at large distances, might have on this result. 

We tested for a density-dependence of the small-scale velocity dispersion, by modeling a linear 
variation of a v with the galaxy density contrast 5 g and determining the coefficient through likeli- 
hood maximization. This significantly improved the VELMOD likelihood, with a best fit relation 
a v = [(100 ± 25) + 35 5 g ] kms -1 . This confirms and strengthens our Paper I result that the galaxy 
velocity field is remarkably cold. Our detection of an increase in a v with density agrees qualitatively 
with the results of Strauss et al. (1998), but our coefficient of 5 g is considerably smaller than their value 
of ~ 50-100 kms -1 . 

We showed that the IRAS-predicted velocity field, with quadrupole, is a good fit to the TF data; the 
correlation function of velocity residuals at f3i = 0.5 is consistent with zero on all scales. Strong velocity 
residual correlations on both small and large scales are seen for f3i < 0.4 and (3i > 0.7, indicating that 
the IRAS-predicted velocity field is not a good fit for these values of f3i. Davis et al. (1996), who adopted 
the Mark III TF zero points, found highly significant discrepancies between the IRAS-predicted and 
Mark Ill-observed velocity fields at all f3i. The VELMOD procedure requires no a priori calibration of 
the TF relation, and with this freedom, the IRAS-predicted velocity field matches the TF data well, 
suggesting that the Davis et al. (1996) discrepancies are tied to uncertainties in the TF calibrations. Our 
claim of agreement between the predicted and observed velocity fields can hold up only if the VELMOD 
TF calibrations ultimately prove correct. 

Indeed, we showed by direct comparison of TF parameters that the VELMOD and Mark III TF 
calibrations (Willick et al. (1997a) differ significantly. The VELMOD TF relations for CF and W91PP 
yield distances ~ 8% shorter than the Mark III TF calibrations, whereas the VELMOD and Mark III TF 
calibrations for A82 and MAT are in good agreement. This has a strong effect on the large-scale bulk 
flow inferred from the data. The VELMOD TF calibrations cannot be brought into closer agreement 
with the Mark III calibration by changing f3j, or by an overall zero point shift in all TF samples. If 
the VELMOD TF relations are correct, then the overall Mark III TF calibration cannot be. Analyses 
based on the published Mark III distances should thus be interpreted with caution. 

The VELMOD TF calibrations are valid, however, only to the degree that the IRAS-predicted peculiar 
velocities are accurate. This will be the case provided that IRAS galaxies trace mass up to linear 
biasing, and linear gravitational instability theory is a good approximation when the galaxy densities 
are smoothed on a 300-500 kms -1 Gaussian scale. Ultimately, the calibration issue must be settled by 
improved observational data. We are carrying out a full-sky TF for this purpose, and will report the 
results of this effort in 1-2 years. 

Our result j3i ~ 0.5 is virtually unchanged from Paper I, ruling out the possibility that cosmic scatter 
and the small volume studied biased our Paper I (3j. Thus, this paper sharpens the discrepancy between 
the VELMOD measurement of f3j and that obtained from the POTIRAS comparison, j3j ~ 0.89 ± 0.12 
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(Sigad et al. 1998). Further underscoring this discrepancy are two analyses that have appeared since 
Paper I, that of Riess et al. (1997) who find fli = 0.4 ± 0.15 using SN la as tracers of the velocity field, 
and that of da Costa et al. (1997) who found /3j = 0.6 ± 0.12 using the SFI TF data set. It may be 
that the differences in the derived values of f3i center on whether the comparison is done at the level 
of the velocities (the v-v comparison, as in this paper, Riess et al., and da Costa et al.) or at the level 
of the densities (the d-d comparison, as in POTIRAS). Future work is needed to determine whether 
these differences can be explained in terms of physical effects, such as a scale-dependent biasing relation 
(e.g., Sigad et al. 1998), or whether they result from TF calibration errors, as discussed above, or other 
methodological factors. The question is an important one because the values of (3i obtained from the v-v 
analyses favor a low-density (Q = 0.2-0.5) universe, while the POTIRAS j3j is suggestive of an Q = 1 
cosmology, if bj £ 1, as suggested by recent analyses of the evolution of rich clusters (Bahcall, Fan, & 
Cen 1997; Fan, Bahcall, & Cen 1997). 

JAW acknowledges the support of NSF grant AST-9617188. MAS acknowledges the support of 
the Alfred P. Sloan Foundation, Research Corporation, and NSF grant AST96-16901. We thank the 
members of the Mark III team, David Burstein, Stephane Courteau, Avishai Dekel, and Sandra Faber, 
for their efforts over the years in putting together the Mark III dataset. We further thank Stephane 
Courteau for discussions concerning selection of the CF sample, and Marc Davis and Tsafrir Kolatt for 
comments on the text. 

REFERENCES 

Aaronson, M., et al. 1982a, ApJS, 50, 241 (A82) 

Aaronson, M., Huchra, J., Mould, J., Schechter, P. L., & Tully, R. B. 1982b, ApJ, 258, 64 

Bahcall, N.A., Fan, X., & Cen, R. 1997, ApJ, 485, 53 

Bertschinger, E., & Dekel, A. 1989, ApJ, 336, L5 

Burstein, D. 1989, privately circulated computer files 

Courteau, S. 1992, Ph.D. Thesis, University of California, Santa Cruz. 

Courteau, S. 1996, ApJS, 103, 363 

Courteau, S. 1997, AJ, 114, 2402 

Courteau, S., Faber, S.M., Dressier, A., & Willick, J.A. 1993, ApJ, 412, L51 

da Costa, L.N., Freudling, W., Wegner, C, Giovanelli, R., Haynes, M.P., & Salzer, J.J. 1996, ApJ, 468, 
L5 

da Costa, L.N., Nusser, A., Freudling, W., Giovanelli, R., Haynes, M.P., Salzer, J.J., & Wegner, G. 
1997, preprint |astro-ph/970T29£ 

Davis, M., Miller, A., & White, S.D.M. 1997, ApJ, 490, 63 

Davis, M., Nusser, A., & Willick, J. A. 1996, ApJ, 473, 22 

Dekel, A. 1994, ARA&A, 32, 371 



30 



Dekel, A. 1997, in Structure Formation in the Universe, eds. A. Dekel and J. Ostriker (Cambridge: 
Cambridge University Press), in press 

Dekel, A., Bertschinger, E., & Faber, S. M. 1990, ApJ, 364, 349 

Dekel, A., Bertschinger, E., Yahil, A., Strauss, M., Davis, M., & Huchra, J. 1993, ApJ, 412, 1 
Fan, X., Bahcall, N.A., & Cen, R. 1997, ApJ, 490, L123 

Fisher, K. B., Huchra, J. P., Strauss, M. A., Davis, M., Yahil, A., & Schlegel, D. 1995, ApJS, 100, 69 

Giovanelli, R., & Haynes, M.P 1985, AJ, 90, 2445 

Giovanelli, R., Haynes, M.P, Myers, S.T., & Roth, J. 1986, AJ, 92, 250 

Giovanelli, R., & Haynes, M.P. 1989, AJ, 97, 633 

Giovanelli, R., Haynes, M., Herter, T., Vogt, N., da Costa, L., Freudling, W., Salzer, J., & Wegner, G. 
1997, AJ, 113, 53 

Han, M.-S., & Mould, J. R. 1992, ApJ, 396, 453 

Hudson, M. J. 1994, MNRAS, 266, 468 

Hudson, M. J., Dekel, A., Courteau, S., Faber, S. M., & Willick, J. A. 1995, MNRAS, 274, 305 
Hudson, M., & Lynden-Bell, D. 1991, MNRAS, 252, 219 

Jackson, J. D. 1976, Classical Electrodynamics, Second Edition (New York: John Wiley) 

Kaiser, N., Efstathiou, G., Ellis, R., Frenk, C, Lawrence, A., Rowan-Robinson, M., & Saunders, W. 
1991, MNRAS, 252, 1 

Kepner, J.V., Summers, F.J., & Strauss, M.A. 1997, NewA, 2, 165 

Mathewson, D. S., Ford, V. L, & Buchhorn, M. 1992, ApJS, 81, 413 (MAT) 

Nusser, A. & Davis, M. 1995, MNRAS, 276, 1391 

Peebles, P. J. E. 1980, Principles of Physical Cosmology (Princeton: Princeton University Press) 

Postman, M. 1995, in Dark Matter, Proceedings of the 5th Maryland Astrophysics Conference, AIP 
Conference Series 336, 371 

Riess, A.G., Davis, M., Baker, J., k Kirshner, R.P 1997, ApJ, 488, LI 

Roth, J. R. 1994, in Cosmic Velocity Fields, eds. F. Bouchet & M. Lachieze-Rey (Gif-sur-Yvette: Edi- 
tions Frontieres), 233 

Schechter, P. L. 1980, AJ, 85, 801 

Schlegel, D. 1995, PhD. Thesis, University of California, Berkeley 

Shaya, E. J., Peebles, P. J. E., & Tully, R. B. 1995, ApJ, 454, 15 

Sigad, Y., Eldar, A., Dekel, A., Strauss, M.A., & Yahil, A. 1998, ApJ, 495, in press 

Strauss, M. A. 1997, in Critical Dialogues in Cosmology, ed. N. Turok (Singapore: World Scientific), 
423 

Strauss, M. A., Davis, M., Yahil, A., & Huchra, J. P. 1992, ApJ, 385, 421 
Strauss, M.A., Ostriker, J. P., & Cen, R. 1998, ApJ, 494, in press 



31 



Strauss, M. A., & Willick, J. A. 1995, Phys. Rep., 261, 271 (SW) 

Tully, R. B. 1988, Nature, 334, 209 

Tully, R. B., & Fisher, J. R. 1977, A&A, 54, 661 (TF) 

Willick, J. A. 1991, PhD. Thesis, University of California, Berkeley 

Willick, J. A. 1994, ApJS, 92, 1 

Willick, J. A., Courteau, S., Faber, S. M., Burstein, D., & Dekel, A. 1995, ApJ, 446, 12 

Willick, J. A., Courteau, S., Faber, S. M., Burstein, D., Dekel, A., & Kolatt, T. 1996, ApJ, 457, 460 

Willick, J. A., Courteau, S., Faber, S. M., Burstein, D., Dekel, A., & Strauss, M. A. 1997a, ApJS, 109, 
333 

Willick, J.A., Strauss, M.A., Dekel, A., & Kolatt, T. 1997b, ApJ, 486, 629 (Paper I) 
Yahil, A., Strauss, M. A., Davis, M., & Huchra, J. P. 1991, ApJ, 372, 380 

A. Appendix: Derivation of Approximate Likelihoods 
The full expressions for the VELMOD likelihoods are given by equations 11 and 12 of Paper I: 



J^°dr r 2 n(r) P(cz\r) S(m } 77, r) expf — ^— ^ M ^) 2 + ^ r ^ 

P(m\r],cz) = ^ — — J ~^\ (Al) 

f£°dr r 2 n{r) P{cz\r) J^dm S (m, i], r) expf — ^l^MM^vM 



f£°drr 2 n(r) $(m — /i(r)) P(cz\r) S{m,ri,r) exp| — - - (^MO)] 

P(j]\m, cz) 



f£°dr r 2 n(r) $(m — fi(r)) P(cz\r) J^dr] S(m,rj,r) exp(— v°(m^(r))] \ (A2) 
where 

, x 1 / \cz - (r + u(r))} 2 \ , K , 

P(czr) = "7=^ exp -I \ 1 >>l , A3 



'V \ V 

S(m,r],r) is the selection function, and /x(r) = 51ogr is the distance modulus. Note the typographical 
error in equation 11 of Paper I; equation ([A1Q is correct. In this Appendix, we derive analytic approxi- 
mations to equations (|A1|) and (|A2|) using the method of steepest descent. We first consider the simple 
case of no selection (S = 1) in A.l, and then consider distance-independent selection functions (A. 2). 
The MAT sample selection function does have a distance dependence; we treat this case in A. 3. In A. 4 
we further refine the approximation and summarize results. 



This piano tables was prepared with the AAS IATgX macros v4.0. 
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A.l. The Case of No Selection 



Equation ( |A3| ) gives the probability that an object at distance r exhibits redshift cz. That probability 
is greatest for r = w, where w is the "crossing point" defined implicitly by cz = w + u(w). Expanding 
about the crossing point gives [cz— (r + u(r))] « — (r — w)(l+u'), where v! is the radial peculiar velocity 
derivative at the crossing point. To the same order of approximation we may write (r — w) ~ w ln(r/w). 
With these approximations, equation (|A3|) becomes: 



P(cz\r) 



1 



exp 



ln(r/w) 
2AT 



(A4) 



where A^ = a v /[w(l + u')]. 

This approximation is valid under certain conditions: First, there must be a unique crossing point 
w. Second, u(r) must be adequately linear within a few times a v of w. Third, w must be sufficiently 
large that the approximation (r — w)/w ~ ln(r/to) is a good one for r within a few times <t„ of w. The 
second and third conditions are satisfied when A v <^1; in practice we found that A v < 0.2 was usually 
sufficient to ensure good accuracy (after the refinements discussed in § A. 4). 

We consider first the forward TF likelihood, equation (|A1|) , in the case of no sample selection (S = 1). 
Substituting equation (Kl) into equation (pO) gives 



P(m\rj, cz) 



/ °° dr r 2 n(r) exp (- 



2A 2 



exp 



\n(r/df 
2A 2 



27T(j TF Jq 00 drr 2 n(r) exp (- 



ln(r/tu) 2 
2A?, 



(A5) 



where Atf = ^V^cttf and d is the forward TF distance (§ 2.3). The integrals in equation [A5] can be 



evaluated analytically if we assume that the density field behaves locally as a power law, 

n(r) = n(w) (—) ■ 
\wj 



(A6) 



In practice, n(r) is not a true power law and the exponent is evaluated as j(w) = [d In n(r)/d In r] 
With this assumption, equation K5| may be written 



P(m\i], cz) 



3 (3+7)3;, 



^ e 



(A7) 



where x = ln(r/w) and y = \n(d/w). The numerator and denominator integrals of equation (A~?) may 
be straightforwardly evaluated to obtain 



P(m\r], cz) 



In 10 



'2nA f 



exp 



2A2 



y-(3 + 7 )A^ 



where 



A, 



A 2 _i_ A2 
^A TF -+- L± v 



1/2 



(A? 



(A9) 



Equation (|A^) has a simple interpretation. When sample selection is neglected, the TF distance d is 
log-normally distributed; the expectation value of In d is In w + (3 + 7) A 2 ,. The fact that ^(lnd) 7^ lnw 
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is due to the Malmquist bias associated with velocity noise; there is both a homogeneous (3A^) and an 
inhomogeneous (jA^) term. Unlike the Malmquist bias in a Method I approach which scales as A FF 
(cf. SW), the bias here is proportional to A^ oc (a v /w) 2 , which is generally much smaller, and which 
decreases with distance. 

The expression for the inverse probability, equation (|A2|) , is complicated by the presence of the 
luminosity function $ in both numerator and denominator. However, like the density field, this function 
varies slowly on the scale relevant to the integration. Consequently, we may treat it too as a power law 
for r near w : 

$(m-//(r)) « $(m-/i(to)) (j^) . (A10) 

Again, we evaluate the power-law exponent according to X(w) = {dln[<&(m — ji(r))]/dhxr} r=w . Once 
this is done, the integrals simplify in the same way as for the forward relation, and we find after similar 
manipulations 

In 10 1 1 f 1 / . , , . 9 \2l 

P(r l \m,cz) = — —f= exp — (|/ illv -(3 + 7 + A)Ay . (All) 

Here y- im = \n(di nv /w), where d inv is the inverse TF distance and e is the inverse TF slope (§ 2.3). The 
fractional inverse TF distance error is given by 



A e i nv 



1/2 



(A12) 



where A^ = (lnl0/5)cr 7? /e. 



Comparison of Eqs. [A8| and |A11| reveals the close analogy between the forward and inverse probability 
expressions when selection is neglected. Such an analogy must indeed hold, for the two forms of the 
TF relation contain the same information. The factor e _1 in equation [All| simply renormalizes the 



probability density to 77-space, while the A reflects the luminosity function dependence of the inverse 
expression. 



A. 2. The role of selection 

In this section, we assume that the sample selection function has no explicit r-dependence, i.e., 
S = S(m,rj). We assume the sample to be selected on a quantity £ with limiting value which is 
linearly related to the TF observables: 

£(m, if) = a\ — bim — c\r\ with scatter erg . (A13) 

The quantities ai, b%, and c\ and <Tg were determined empirically for the Mark III samples by Willick et 
al. (1995, 1996). Willick (1994) shows that: 

S(m,r)) = ~ [1 + erf (At(m,r))} , (A14) 



where 



^ (m ,,) 3 «^Z«£. (A15) 
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We define a TF-predicted apparent magnitude m r = M(rj) + 5 log r. Then, using the identities derived 
by Willick (1994), the forward likelihood becomes: 



P(m\r], cz) 



[1 + erf (y4 ? (m, 77))] 



J r 2 n(r)P(cz\r) exp ^— ^ rfr 



27T(T TF / r 2 n(r)P(cz\r) 1 + erf (A^(m r ,r])/^1 + /? 2 ) 



(A16) 



where /3 = &i<7tf/°"£- 

The integral over m has caused the A$ term in the deonominator to acquire an r-dependence, although 
it did not start out with one. This complication makes it inconvenient to follow our previous procedure 
exactly. Instead, we treat this term as constant across the effective range of integration, and take it 
outside the integral; this is correct to the same order of approximation. This leaves us with a ratio of 
integrals we have already evaluated. We then require that the resultant probability density P(m\r), cz) 
be properly normalized, yielding: 



P(m\r), cz) = 



In 10 1 + erf [A^(m, r})] 



erf 



^-p{-^|0<-(3 + ^) 2 }, (A17) 



where 



and 



m (r],w) = M(rj) + 5 log w + [3 + 7] K 



(3 = 



where a e = 5A e / In 10 . 



(A18) 



(A19) 



The effect of selection appears purely outside the exponent now. Indeed, the role of selection is very 
similar to what it was in pure Method II (Willick 1994), with a slightly different evaluation of A$ and 
j3 in the denominator. 

The corresponding expression for the inverse relation follows directly, given the analogy we drew 
between the two expressions in the previous subsection: 



P(r)\m, cz) = 



In 10 1 + erf [A^{m, r})] 1 



1 + erf 



2tt A, 



eXP {"2A^( y_(3 + 7 + A)A ') 2 } ' 



where 



and 



r] (m, w) = rf [m — 



5logw+ hTTo (3 + 7 + A)A ' 



P = where a v e = Jal + (5e/ In 10A„) 



(A20) 



(A21) 



(A22) 



Note the different definition of (5 in the inverse and forward cases. In particular, if selection is r\- 
independent ic\ = 0), the terms involving the error functions cancel, and P(r]\m, cz) reduces to the no 
selection expected. 
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A. 3. Treating an explicitly distance-dependent selection function 

If the selection function S has an explicit distance dependence, things get a bit more complicated. In 
Willick et al. (1996), the data for all the Mark III samples was fit to the form, 

£ = £ (m, 77, r) = di — b\m — C\r\ — d\ log r ; (A23) 

only MAT had a significantly non-zero value of d\. However, c\ = for MAT; selection for MAT has no 
explicit ^-dependence and we take this into account in what follows. Corresponding to equation (|A23 ) 
is an r-dependent A% parameter, 



A$(m, r) 



£(m, r) - & 



(A24) 



and thus the selection function S(m,r) = [1 + erf (A^(m,r)] /2. 

The main effect of distance-dependent selection is to introduce a new power-law exponent, a, into 
our earlier expressions, where 



a(m, w) = 



dlnS 



dliar 



di 



7T In 10 a i [1 + erf (A^)] 



(A25) 



For the inverse relation, the addition of a is all that is required to correct our expressions. Specifically, 



P(j]\m, cz) 
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1 



1 



exp 



e,vnv 



2A 2 ■ 

e,mv 



Vinv - (3 + 7 + ^ + Oi)A 



(A26) 



There are no selection functions out in front because for MAT, selection is r\- independent. 

For the forward relation, the fact that a depends on m ruins the pure Gaussianity of the exponent. 
Using the same approach as we did to derive equation ( |A17|) , we assume that a varies slowly with m, 
take Ae out of the denominator integral, and normalize after the fact. After some algebra, one finds 



P(m\r], cz) 



In 10 1 + erf (A^(m,w)) 



1 + erf 
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2,71 A' 



exp 



^(l/-(3 + 7 + «(m, W ))A 2 u )' 



2AI 



(A27) 



In equation |A27|, the individual terms have the following definitions: 



m = m + a(m , w) 



——A 2 
In 10 u ' 



(A28) 



1 - 



5 . 9 da 

A 2 — 

In 10 v dm 



(5 = 



5V2b 1 a(m,w)A 2 v ( 
lnlOfT, \ i 



where a' 



sa; 

In 10 



a(m, w) In 10 
V2d~i 



n -1 



(A29) 
(A30) 
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A. 4. Final Refinement and Summary 

Our original approximation to P(cz\r), equation ( |A4[ ), was correct to first order in a v . This leads 
to systematic inaccuracies in two regimes: small distances (,$,2000 kins" 1 ), where the approximation 
(r — w)/w ~ \n(r/w) loses accuracy, and in regions of velocity field curvature, when u"a v is comparable 
to u' . We extend its regime of validity by making second-order corrections for these effects. 

To second order in x = \n{r/w) we find (r — w)/w = x + x 2 /2. Using this and a second order Taylor 
expansion of u(r) about w, and retaining only terms of order (a v /w) 2 in the exponent, we find after 
some algebra 



where 



exp 



2^ 2 



(cz — [r + u(r)]Y 



u"w 



1 + u' 



x f(r) 



(A31) 



(A32) 



where u' and u" are evaluated at w. The term e x ^ 2Av in equation ( |A31[ ) is just our original approxima- 



tion for P(cz\r), equation ( |A4j) , and f(r) is the second-order correction. It is non-Gaussian in ln(r/w) 
and thus cannot be analytically integrated as before. 

We thus treat it as we have other slowly-varying terms: we approximate it as a power law in the 
vicinity of the crossing point. However, because of its cubic nature, the local logarithmic derivative 
is identically zero. We thus proceed heuristically by calculating the power-law exponent as a finite 
difference over an interval of lnr of ±gA v , where g is of order unity: 



v 



In f(x = gA v ) - In f(x = -gA v ) 
2gA v 



9 



(A33) 



We calibrated the appropriate value of g by varying it until we maximized agreement between the exact 
and approximate likelihoods. This happened at g — 1.5, and thus the correct exponent is 



-1.1 (1+e) 



(A34) 



This leads to the final forms of the analytic approximation to the VELMOD likelihoods. For the 
forward relation, P(m\r], cz) is given by equation ( A17 ) for A82, W91PP, and CF (the samples for which 
the selection function has no explicit distance dependence) and by equation (|A27j ) for MAT. For the 
inverse relation, P(r]\m, cz) is given by equation (|A20|) for A82, W91PP, and CF and by equation flA"2l| 
for MAT. However, in all of these equations, the quantity 3 + 7 is replaced by 3 + 7 + u, where v is 
given by equation ( |A34j) , and e is given by equation (|A32j ). Note that the definition of v is such that 
the homogeneous Malmquist bias term is reduced from 3A^ to 1.9 A 2 . This is a significant effect for 
distances ,$,2000 kms" 1 , and thus the refinement discussed here is crucial for extending the regime of 
validity of the approximation to small distances. 
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